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ABSTRACT 

We have searched for unresolved X-ray sources in the vicinity of two rich clusters of galaxies: Abell 1995 
(A1995) and MS 0451.6-0305 (MS0451), using the Chandra X-ray observatory. We detected significantly more 
unresolved sources around A 1995 than expected based on the number of X-ray sources to the same flux limit 
detected in deep Chandra observations of blank fields. Previous studies have also found excess X-ray sources in 
the vicinity of several nearby clusters of galaxies using ROSAT , and recently in more distant (z rs 0.5) clusters 
(RXJ0030 and 3C295) using Chandra . In contrast, we detect only 14 unresolved X-ray sources near MS0451, 
which is consistent with the number expected from a cluster-free background. We determine the luminosity 
functions of the extra sources under the assumption that they are at the distance of their respective clusters. The 
characteristic luminosity of the extra sources around A 1995 must be an order of magnitude fainter than that of the 
extra sources around RXJ0030 and 3C295. The apparent lack of extra sources around MS0451 is consistent with 
its greater distance and the same characteristic luminosity as the A1995 sources. Hardness ratios suggest that, on 
average, the extra sources in A 1995 may have harder spectra than those of RXJ0030 and 3C295. These results 
indicate that different classes of objects may dominate in different clusters, perhaps depending on the formation 
history and/or dynamical state of the accompanying cluster. 

Subject headings: galaxies: clusters: individual (Abell 1995, MS 0451.6-0305) — X-rays: galaxies: clusters 



1. INTRODUCTION 

Evidence has accumulated recently that there are more X-ray 
point sources in the direction of clusters of galaxies than to- 
ward cluster- free regions of the sky. Henry & Briel (1991), us- 
ing ROSAT PSPC observations, found just about twice as many 
unresolved sources around Abell 2256 (at z = 0.06, Struble & 
Rood 1991) as expected from blank field (no clusters) back- 
ground observations. The luminosity of these sources, assum- 
ing they are at the redshift of the cluster was found to be about 
10 42 erg s -1 or greater (in 0.5 - 2 keV). These sources have high 
X-ray to optical flux ratios. Some of the sources in A2256 were 
identified as cluster member galaxies. Henry and Briel also dis- 
cuss the possibility that the emission from these sources is due 
to hot gas in galaxies not removed by ram pressure or evapora- 
tion, or due to shocks in gas from merging. 

Similarly, Lazzati et al. (1998), analyzing ROSAT PSPC im- 
ages, found an excess number of unresolved X-ray sources in 
the fields of two nearby clusters: A194 and A1367 (z = 0.018 
and 0.022). The spectra of the sources were consistent with 
thermal bremsstrahlung with T < 2 keV. Lazzati et al. also 
found evidence for association of some of these sources with 
cluster member galaxies, implying luminosities between 0.6 
and 6.6 xlO 41 erg s" 1 in the 0.5 - 2 keV band. X-ray emis- 
sion from hot gas associated with cluster member galaxies had 
been reported earlier, also based on ROSAT PSPC observations 
(Grebenev et al. 1995; Bechtold et al. 1983). 

Indirect evidence has also been presented for the existence 
of an excess population of unresolved X-ray sources associated 
with Abell clusters. Soltan & Fabricant (1990) using Imaging 
Proportional Counter data from the Einstein Observatory found 
excess fluctuations in nearby galaxy clusters which could be 
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explained by assuming the presence of low luminosity sources 
(«4x 10 41 erg s -1 ) in clusters with extent less than 1'. They 
discuss the possibility that the emission from these sources is 
due to low luminosity AGNs, or to hot gas in member galax- 
ies. Soltan etal. (1996) found a correlation between the surface 
brightness of the X-ray background and Abell clusters on scales 
of a degree, which is much larger than the X-ray emission from 
the intracluster gas. The characteristic length was found to be 
about 10 h~ l Mpc in radius (where H = /il00 Km s"' Mpc" 1 ), 
i.e. extra X-ray emission was found around Abell clusters out 
to about 15 Mpc (h = 0.65). Soltan et al. could not explain the 
extra X-ray emission based on known sources, or random fluc- 
tuations in their number density. They estimated the required 
number of excess sources to be about 50% above the expected 
number of background sources. 

Most recently Cappi et al. (2001), using Chandra ACIS (Ad- 
vanced CCD Imaging Spectrometer) observations, found twice 
as many unresolved X-ray sources in the images of two distant 
clusters of galaxies, 3C295 (z = 0.46 Dressier & Gunn 1992), 
and RX J003033.2+261819 (RXJ0030; z = 0.5 Vikhlinin et al. 
1998), as expected from a cluster- free background to their flux 
limits (Giacconi et al. 2001; Mushotzky et al. 2000). 

Our main goal in this letter is to present new results on the 
flux and number distributions of unresolved X-ray sources based 
on our Chandra ACIS observations of two rich clusters of galax- 
ies: A1995 and MS0451. We also address briefly the nature of 
the excess sources. 

2. data processing and analysis 

A 1995 is a rich cluster at z = 0.32 with an intracluster gas 
temperature of T x = 7.6 keV (Patel et al. 2000). A1995 was 
observed with the Chandra ACIS in May and July 2000 for 35 
k sec and 12 k sec. The aim point was on the back-illuminated 
chip S3 of ACIS-S. The data were taken in full frame mode with 
a readout time of 3.2 sec. We used standard Chandra software 
tools to clean the data of time intervals with high background 
and/or bad aspect, remove bad or flickering pixels, and correct 
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Table 1 

Unresolved X-ray sources detected in A1995. 



X-RAY SOURCE 


COUNTS 


COUNTS 


R" 




0.5-2 KEV 


2-10 KEV 


Mag 


CXOU J45308.70+580313.8 


585±25 


179 ±14 


18.8 


CXOUJ45305. 82+580309.1 


342±19 


79.1 ±9.6 


19.2 


CXOU J45307. 10+580205. 8 


306±18 


74.2 ±9.2 


20.8 


CXOU J45327.65+580339.5 


166±14 


39.7 ±7.6 


20.9 


CXOU J45305.45+580033.9* 


69.5±8.5 


<1() 


10.2 


CXOU J45246.31+580059.7 


65.4±8.2 


14.5 ±4.0 


20.9 


CXOU J453 17.46+580003.0 


56.3±8.2 


20.9 ±5.9 


20.3 


CXOU J45233. 22+580559.2 


48.8±7.1 


17.6 ±4.4 


19.5 


CXOU J45229. 56+580418. 2 


37.5±6.2 


20.4 ±4.7 


>22 


CXOU J45230.74+580448.5 C 


34.5±6.0 


<10 


16.0 


CXOU J45233. 56+580456. 3 


28.4±5.4 


23.0 ±5.0 


22.0 


CXOU J45324. 68+580318.4 


26.7±5.7 


<1() 


>22 


CXOU J45316. 75+575928. 6 


26.6±6.1 


<10 


>22 


CXOU J45248. 65+580255. 5* 


24.3±5.3 


<10 


13.2 


CXOU J45301. 83+580005. 7 


22.1±5.0 


11.4 ±4.0 


>22 


CXOU J45255. 19+580056.0 


20.7±4.7 


<10 


>22 


CXOU J45 244. 00+5 80203. 3 


20.6±4.7 


14.2 ±4.0 


>22 


CXOU J453 15.41+580448.5 


20.5±5.6 


<1() 


21.0 


CXOU J45253.64+580020.3 


20.1±4.6 


<10 


>22 


CXOU J45319.1 1+580134.4 


17.2±4.8 


<1() 


>22 


CXOU J45313.22+580126.9 


16.5±4.6 


<10 


>22 


CXOUJ45315.55+580117.4 


16.3±4.9 


<1() 


20.3 


CXOU J45242.50+580159.1 


15.8±4.1 


<10 


>22 


CXOU J45228.20+575954.5 


15.0±4.0 


7.2 ±3.0 


21.1 


CXOU J4523 1.03+580010.9 


14.5±3.9 


5.1 ±2.4 


21.1 


CXOU J45235.73+580656.1 


12.8±3.7 


<1() 


>22 


CXOU J45245.50+580519.5 


12.2±3.6 


<1() 


>22 


CXOU J45234.49+575904.0 


11.9±3.6 


14.2 ±4.2 


>22 


CXOU J4525 1.93+580046.2 


10.1±3.3 


12.3 ±3.7 


20.9 


CXOU J45322.04+575858.7 


<10 


25.3 ±6.4 


>22 



" R band magnitude of optical counterparts/limit if not detected 

* GSC2 object (http : //www-gsss . stsci . edu/gsc/gsc2) 

c IRAS Galaxy (F1451 1+5816) 



for event gains for a focal plane temperature of -120 C. The 
final effective exposure time of the merged data after cleaning 
was 54.5 k sec (for details see Joy et al. 2002). MS0451 is 
a rich cluster at z = 0.55, with an intracluster gas temperature 
of T x = 10.9 ± 1.2 keV (Donahue et al. 1999). MS0451 was 
observed with Chandra in October 2000 for 45 k sec, also on 
the back-illuminated chip S3. After similar data processing, the 
cleaned data set has a 44.75 k sec exposure time (for details see 
Donahue et al. 2002). 

We analyzed data from within < 5' of the nominal center of 
the clusters, close to the optical axis, where the point spread 
function (PSF) is not degraded significantly, and the change in 
effective area is negligible. We made images in the 0.5 - 2 keV 
and 2-10 keV energy bands, where the calibration is most re- 
liable. We used the wavedetect package of the Chandra 
Interactive Analysis of Observations setting the probability of 
fake detection to lxlO" 6 (corresponding to 4.7 a, for details, 
see Cappi et al. 2001). We kept sources which had a signal-to 
noise ratio of 3 (using a local background estimate) correspond- 
ing to the detection of about 10 counts within the source area. 
On average, we would expect one spurious source in an area of 
the sky four times that of the S3 chip. These are similar to the 
conditions used by Cappi et al. (2001). The background clus- 
ter emission is very smooth, and, according to our Monte Carlo 
simulations, the probability of the cluster emission resulting in 
false detection of a point source is negligible. We used a power 
law with a photon index of T = 1 .4 when converting count rates 



Table 2 

Unresolved X-ray sources detected in MS045 1 . 



X-RAY SOURCE 






COUNTS 


COUNTS 








0.5-2 KEV 


2-10 KEV 


CXOUJ45419.63- 


-30420 


.5 a ' b 


746 ±28 


150 ±13 


CXOU J45356.32- 


-25837 


.7° 


396 ±20 


121 ±11 


CXOU J45422.59- 


-30035 


.2 


87.9 ± 9.5 


17.8 ±4.4 


LAUU J4J424. / j- 




o 
.IS 


oy.u ± o.D 


n 1 J-A n 

JJ.L ±O.U 


CXOU J45426.07- 


-30013 


.2 


57.1 ±7.7 


12.1 ±3.7 


CXOUJ45412.81- 


-30047 


.7 


45.2 ± 8.1 


<10 


CXOUJ45419.20- 


-30521 


.2 


28.4 ± 6.0 


<10 


CXOU J45408.57- 


-30521 


.2 


27.4 ± 6.0 


<1() 


CXOU J454 10.88- 


-30125 


.2 


21.4 ±5.7 


<10 


CXOU J45355.65- 


-30409 


.6" 


19.7 ±4.9 


<10 


CXOU J45406.70- 


-30412 


.3°'* 


18.0 ±4.6 


<10 


CXOUJ45421.95- 


-25816 


.2 


17.2 ±4.2 


43.9 ±6.8 


CXOUJ45404.19- 


-30403 


.7 


12.0 ± 3.9 


<10 


CXOUJ45421.39- 


-30132 


.4 


10.9 ± 3.5 


<10 


CXOU J45356.73- 


-30226 


.1 


<10 


23.1 ±5.6 



1 GSC2 object 

2 2MASS source (http : / /www . ipac . caltech . edu/2mass) 



to fluxes. We did not correct for vignetting since in this soft 
band (0.5 - 2 keV) and off axis angles of < 5' the correction 
is negligible. Note, that in this soft band the fluxes have only a 
weak dependence on the power law slope and absorbing column 
density, which are A/ H = 3.5 x 10 20 cm -2 and 3.9 x 10 20 cm -2 for 
A 1995 and MS0451. 

3. RESULTS 

We detected 29 and 14 unresolved X-ray sources in the fields 
of A 1995 and MS0451. The sources show no correlation with 
the spatial distribution of the cluster emission. Source details 
are given in Table 1 and 2 and the logAZ-logS curves are plot- 
ted in Figure 1. In Figure 1 we also show the log N- log S curves 
for the sources near RXJ0030 and 3C295 (Cappi et al. 2001). 
We quote source densities in terms of number per ACIS chip 
(8' x 8') in order to show the actual numbers of unresolved 
sources as detected. The expected logA/- log S curves from 
cluster-free background (dashed lines) is taken from Mushotzky 
et al. (2000) which predicts a slightly higher number density 
of background sources than Rosati et al. (2002), Campana et 
al. (2001) and Giacconi et al. (2001), and slightly less than 
Brandt et al. (2001). We choose the background predicted by 
Mushotzky et al. (2000) since it provides the best fit to the high 
flux end of our logA/-log S curve. We detect sources with fluxes 
brighter than 6 x 10" 16 erg s" 1 cm -2 and 8 x 10" 16 ergs" 1 cm -2 
in A 1995 and MS0451. We obtained R band magnitudes for 
unresolved X-ray sources in A 1995 (see details in Patel et al. 
2000). The 0.5-2 keV X-ray to R band optical flux ratios of 
unresolved sources in the field of A 1995 (cf. Table 1) seem 
to have a similar distribution to those of background sources 
(Mushotzky et al. 2000). 

The logA/-logS curve of the unresolved sources in the A 1995 
field is steeper than that of the background between 1 and 3 x 
10~ 15 erg s" 1 cm" 2 , which indicates a build up of extra sources 
(excess number of sources relative to the background) with fluxes 
in this interval. There is an indication that for fluxes less than 
10~ 15 erg s -1 cm" 2 the slope of the logAZ-logS curve of A1995 
is close to that of the background, suggesting that there are no 
extra sources with flux below this value. The logA/-log S curves 
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of unresolved sources in 3C295 and RXJ0030 show similar 
signs of a cut off at low fluxes in the distribution of the excess 
sources (cf. Figure 1). Overall we find 29 unresolved sources 
in the A 1995 field, as opposed to the expected number for a 
cluster free background of about 17. 

In the field of our more distant cluster, MS0451, we find 14 
unresolved sources, which is within la from the number ex- 
pected based on a cluster free background. Since MS045 1 is 
at about the same redshift as RXJ0030 and 3C295, our expo- 
sure time is longer than those of RXJ0030 and 3C295, and we 
found no extra sources in the MS045 1 field, we have the sim- 
ple and potentially important finding that: not all clusters have 
extra unresolved sources associated with them at the flux limits 
of these Chandra observations. 




S [CGS] 



FIG. 1.— log/V(> S)-logS per ACIS chip (8' X 8') in the 0.5-2.0 keV 
band for unresolved X-ray sources in the fields of Abell 1995 and MS0451 
(squares and triangles, our results), and RXJ0030 and 3C295 (stars and dia- 
monds, Cappi et al 2001). The expected background is shown with dashed 
line. The expected background plus cluster sources, using their derived lumi- 
nosity functions (see text), are shown using solid lines. On the left Y axis we 
show the corresponding numbers per square degree for comparison. 

In this letter we estimate the luminosity functions for the ex- 
tra sources in A1995, RXJ0030 and 3C295 using their logN- 
logS curves, assuming the excess sources are at the redshifts 
of their respective clusters. As customary, we used a lumi- 
nosity function of the form $(L) oc (L/L^)~ a with a lower cut 
off at a characteristic luminosity, L* (i.e. there are no sources 
with L < L*). The fits, which are a good description of the 
logiV — log S distribution, are shown using solid lines in Fig- 
ure 1. We found that the same slope, a = 3.1, could describe 
the unresolved sources in all three clusters. However, the char- 
acteristic luminosities were found to be L*= 0.5, 4.0, and 4.8 
x 10 42 erg s" 1 in A1995, RXJ0030 and 3C295. The total num- 
ber of extra sources based on these fits are: 10, 12, and 7 (with 
about ±3 statistical error) for A1995, RXJ0030 and 3C295. 
Due to the cluster emission, our ability to detect faint sources 
decreases toward the center (about 20 counts per detect cell for 
A1995). This effect is a fraction of the Poisson fluctuations 
(it corresponds to missing at most a single additional source) 
and so we ignore it. We verified, by means of Monte Carlo 
simulations, that the observations can be drawn from the as- 
sumed background distribution plus a distribution based on the 
derived luminosity functions of the extra cluster sources, and 
that the observations can not be explained using the same lu- 
minosity function for all clusters. Although the uncertainties 
in the luminosity function parameters are quite high: ±0.5 in 
the slope, and 30%-40% in the lower cut-off, the character- 
istic luminosity of unresolved sources associated with A 1995 



are significantly (about one order of magnitude) less than those 
associated with RXJ0030 and 3C295. Furthermore, when the 
luminosity function of unresolved sources near A 1995 scaled 
to the greater distance of MS0451 the resulting logN-logS is 
fully consistent with the observed one (see Figure 1, solid curve 
associated with MS0451). On the other hand, it is not consis- 
tent with the scaled luminosity function of either RXJ0030 or 
3C295 (which are roughly the same distance as MS0451). We 
expect the effect of spatial variations in the PSF, and effective 
area, on the completeness function to be small, with an over- 
all effect on the derived luminosity function much less than the 
errors quoted above. 




0.0001 0.0010 0.0100 0.1000 

H+S[CTS/SEC] 

FIG. 2. — Hardness ratios of unresolved X-ray sources as a function of 0.5- 
10 keV count rate in the A1995, MS0451, RXJ0300, and 3C295 fields. The 
symbols are the same as in Figure 1 . 

4. DISCUSSION 

The derived luminosity functions enable us to estimate the 
contribution of the unresolved sources to the overall X-ray emis- 
sion around clusters of galaxies. The derived surface brightness 
of unresolved X-ray sources in the 8' x 8' fields near RXJ0030 
and 3C295 is comparable to the surface brightness of the X- 
ray background (2.6 x 10~ 8 erg s -1 cm" 2 str _1 in the 0.5-2 keV 
band). This enhancement is about 100 times larger than the cen- 
tral enhancement from large scale emission (10 /T 1 Mpc) found 
by Soltan et al. (1996). Cappi et al.'s results (cf. their Figure 
4) show that these unresolved sources do not extend beyond 
the clusters much more than about 4', as opposed to the large 
scale component of Soltan et al., which extends out to about 30' 
when scaled to the redshifts of RXJ0030 and 3C295 (z ~ 0.5). 
Therefore it is likely that this component contributes only to the 
compact component found by Soltan et al. 

There are a number of possibilities for the nature of these 
extra sources: cosmic variance, star-burst galaxies, a result of 
gravitational lensing of background objects, or an enhanced 
number density of AGNs/QSOs. It is unlikely that the extra 
sources are due to cosmic fluctuations, which is only at the level 
of about 20%-30%, significantly below the measured factor of 
two (Cappi et al. 2001). It could be possible, however, that 
the excess of unresolved sources is due to projection effects, 
with differences arising from whether we are viewing along or 
perpendicular to a filament of the cosmic web. Star-burst galax- 
ies are also unlikely to be the sources Cappi et al. found since 
their X-ray luminosities are about 10-100 times too faint. How- 
ever, recent Chandra deep surveys did find some exceptional 
X-ray bright galaxies at similar redshifts. Gravitational lens- 
ing can increase the number of unresolved sources, but only if 
the \ogN-\ogS slope is steep enough (> 0.4, Croom & Shanks 
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1999; Mellier 1999). Two opposite effects are competing in de- 
termining the number of observed sources: lensing magnifies 
the flux, but it also reduces the field of view behind the gravi- 
tational lens. Lensing would need a significantly higher slope 
in log N — log S to explain the large number of extra sources in 
the field. Refregier & Loeb (1997) predict an average reduc- 
tion of the surface density of faint sources at fluxes less than 
10" 15 erg s" 1 cm" 2 . Cappi et al. conclude that, as far as spec- 
tra and luminosities are concerned, the unresolved sources near 
RXJ0030 and 3C295 could be AGNs/QSOs associated with the 
respective clusters. 

In Figure 2 we show the hardness ratios (HRs), (H-S)/(H+ 
S), where S and H are X-ray fluxes in the 0.5-2 keV (soft) and in 
the 2-10 keV (hard) bands, of unresolved sources in the fields 
of A 1995 (squares), MS045 (triangles) as a function of H + S 
(our results). As a comparison, we also plot the hardness ratios 
of RXJ0300 (stars) and 3C295 (diamonds) from Cappi et al. 
(2001). Points with one sided error bars represent sources not 
detected either in the hard or in the soft band. The average 
HR of unresolved sources detected in both soft and hard bands 
near RXJ0300 and 3C295 are w -0.5, while the average HR of 
unresolved sources near A1995 is slightly harder, about -0.25. 
These sources near A 1995 also have lower fluxes than sources 
near the other two clusters, as previously noted (see section 3). 

Comparing our Figure 2 to Figure 3 of Rosati et al. (2002), 
which shows the HRs of sources of different types as a func- 
tion of their luminosities, we conclude that unresolved sources 
in RXJ0300 and 3C295 with H+S w 0.004 cts/sec, correspond- 
ing to luminosities of about 10 44 erg s" 1 at the distance of the 
clusters, would be compatible to the HRs of Type I AGNs (as 
noted by Cappi et al. 2001). While the faint unresolved sources 
(L* = 5 x 10 41 erg s" 1 ) in the field of A 1995 are concentrated 
around H + S « 0.0006 cts/sec, HR « -0.25, which falls be- 
tween normal and star-burst galaxies. 

Recent results show that the angular correlation function of 
X-ray selected AGNs is similar to that of nearby galaxies sug- 
gesting that AGNs sample the mass density the same way as 
galaxies sample (Akylas, Georgantopoulos & Plionis 2000), 
in contrast to optically selected AGNs, which are found to be 
more frequent in field galaxies (5%) than in galaxies near clus- 
ters (1%, Dressier Thompson and Shectman 1985; Osterbrock 
1960). Therefore we would expect more X-ray selected AGNs 
in clusters. Since the clustering length of X-ray selected AGNs 
and nearby galaxies is the same within errors (rs 7 h~ l Mpc, 
Basilakos 2001; Akylas, Georgantopoulos & Plionis 2000; Pee- 
bles 1993), we would expect the ratio of total number of galax- 
ies to the number of X-ray selected AGNs, N ga i /Nagn, to be 
~ (n a i)/(nAGN) (where (n ga i) and («agw) are the average num- 
ber densities of galaxies and AGNs in the cluster). However, 
this effect could only account for about 20% of the excess, 
much less than the a factor of two, which has been found by 
Cappi et al (2001) and this work. 

Our results, that the characteristic luminosities of extra sources 
are about one order of magnitude different in A1995 vs. 3C295 
and RXJ0030, and that there seems to be a difference between 
their HRs argues against cosmic variance and projection effects. 
It suggests instead that different class of objects might dominate 
in different clusters perhaps depending on the formation history 
and/or the dynamical state of the cluster. 

Perhaps the unresolved sources in A 1995 belong to a class 
of starburst galaxies, a result of enhanced star formation due to 
interactions between infalling groups of galaxies and the intra- 



cluster gas. This enhanced star formation would lead to an ex- 
cess of blue galaxies around these areas similar to the Butcher- 
Oemler effect (Butcher & Oemler 1978). A search for a cor- 
relation between galaxy color changes around X-ray sources 
compared to other areas in the cluster could be used to check 
this possibility. 

At present, the exact nature of these objects is not known. 
Due to limited photon statistic, their spectra could not be de- 
termined individually, both low temperature (< 2 keV) ther- 
mal bremsstrahlung and power low spectra can be fitted to their 
stacked spectra. Revealing the physical properties of these ob- 
jects would help us to improve our understanding of structure 
formation, specifically the origin and evolution of the intra- 
cluster gas, and the effect of merging. Identification of these 
sources would also help to asses the contamination these sources 
cause in the interpretation of cluster emission as thermal brems- 
strahlung from intra-cluster gas. This contamination would re- 
sult an overestimation of the normalization of the X-ray flux 
from the cluster and would lead to a systematic error in the 
determination of the Hubble constant using SZ effect and ther- 
mal bremsstrahlung (see for example: Molnar, Birkinshaw and 
Mushotzky 2002). Follow up observations of the individual 
sources are necessary to solve this mystery. 

This work was partly supported by NASA LTS A Grant NAG5- 
3432 and Chandra grant GO0-1049C. The optical analysis was 
done by Tom Brink. We acknowledge useful discussions with 
Mark Birkinshaw, Renyue Cen, Richard Mushotzky and Craig 
Sarazin. 
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